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. Abstract 

'""i ' By constructing a coupling with unbounded time-dependent drift, dimension- 

. free Harnack inequalities are established for a large class of stochastic differential 

Q^ I equations with multiplicative noise. These inequalities are applied to the study of 

■ heat kernel upper bound and contractivity properties of the semigroup. The main 

■ results are also extended to reflecting diffusion processes on Riemannian manifolds 
^ I with non-convex boundary. 

1 Introduction 

Consider the following SDE on W^: 

(1.1) dX, = a(t, Xt)dBt + b{t, Xt)dt, 

where Bt is the ci-dimensional Brownian motion on a complete filtered probability space 
(fi,{^t}t>o,P), and 
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a : [0, cx)) X M'^ X ^ M'^ ® R'^, b : [0, oo) x R"^ x ^ M"' 

are progressively measurable and continuous in the second variable. Throughout the 
paper we assume that for any Xq G M*^ the equation (11. ip has a unique strong solution 
which is non-explosive and continuous in t. 

Let Xf be the solution to (11. ip for Xo = x. We aim to establish the Harnack inequality 
for the operator Pt. 

PJix) := E/(Xf), t > 0,x G M^/ G ^+(M'), 

where J3§^{M.'^) is the class of all bounded non-negative measurable functions on M'^. To 
this end, we shall make use of the following assumptions. 

(Al) There exists an increasing function K : [0, oo) — ?■ R such that almost surely 

Wait, x) - ait, y)\\ls + 2(K^, x) - bit, y)^x - y) < Kt\x - y\\ x,y e R^ t > 0. 

(y42) There exists a decreasing function A : [0, oo) — (0, oo) such that almost surely 

ait, x) V(t, x) > xeW^,t>Q. 

(y43) There exists an increasing function 5 : [0, oo) — (0, oo) such that almost surely 

\iait,x)-ait,y))ix-y)\<5t\x-y\, x,y eW,t>Q. 
(y44) For n > 1 there exists a constant c„ > such that almost surely 

\\ait,x) - ait,y)\\Hs + W^^) - < Cn|x - y\, \x\, \y\,t < n. 

It is well known that iAl) ensures the uniqueness of the solution to (11. ip while (^4) 
implies the existence and the uniqueness of the strong solution (see e.g. [12] and references 
within for weaker conditions). On the other hand, if b and a depend only on the variable 
X G R"', then their continuity in x implies the existence of weak solutions (see [T^ Theorem 
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2.3]), so that by the Yamada-Watanabe principle [28], the uniqueness ensured by (Al) 
imphes the existence and uniqueness of the strong solution. 

Note that if a(t, x) and b{t, x) are deterministic and independent of t, then the solution 
is a time-homogeneous Markov process generated by 



L := ^ ^ ttijdidj + ^ bidi, 

where a := aa*. If further more a and b are smooth, we may consider the Bakry- Emery 
curvature condition [5]: 

(1.2) r2(/,/)>-i^r(/,/), feC^iR') 

for some constant G M, where 



1 

r(/,^?) := 2 E <^ijW){d,g), f,g G C'iR"), 

T,{f, /) := hvif, /) - r(/, L/), / G C^^IM-^). 

According to [221 Lemma 2.2] and [231 Theorem 1.2], the curvature condition (11. 2p is 
equivalent to the dimension- free Harnack inequality 



{Ptf{x)r<{PtF{y))exp 
where 



2(p- 1)(1 -e-2^*) 



, t>0,p>l,/G^+(M'^),x,2/G 



p,(x,|/):=sup{|/(x)-/(y)|: / G ^^(M'^), r(/, /) < l}, x,y eW. 

This type of inequality has been extended and applied to the study of heat kernel (or 
transition probability) and contractivity properties for diffusion semigroups, see [H [THl S] 
for diffusions on manifolds with possibly unbounded below curvature, [2^ [T6] for stochas- 
tic generalized porous media and fast diffusion equations, and [2 [3l El HZl [iSl [191 ttB 129] 
for the study of some other SPDEs with additive noise. 
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If a depends on x, however, it is normally very hard to verify the curvature condition 
(11 ■2p . which depends on second order derivatives of a~^, the inverse matrix of a. This is 
the main reason why existing results on the dimension-free Harnack inequality for SPDEs 
are only proved for the additive noise case. 

In this paper we shall use the coupling argument developed in j^, which will allow us 
to establish Harnack inequalities for cr(t, x) depending on x. This method has also been 
applied to the study of SPDEs in the above mentioned references. To see the difficulty in 
the study for o"(t, x) depending on x, let us briefly recall the main idea of this argument. 

To explain the main idea of the coupling, we first consider the easy case where a and 
h are independent of the second variable. For x ^ y and T > 0, let Xt solve (11. ip with 
Xq = X and Yt solve 

dy, = a{t)dB, + h{t)dt + t^M(^i_^dt, Fo = y- 

Then Yt is well defined up to the coupling time 



r:=inf{t>0: Xt = Yt}. 



Let Xt = Yt for t > r. We have 

d\Xt~Yt\ = -^-^jr^dt, t<T. 

This implies r = T and hence, Xt = Yt- On the other hand, by the Girsanov theorem 
we have 



Prfiy) = ElRfiYr)] 



for 



R := exp 
Therefore, 



x-y\ {a{t)-\Xt-Yt),dBt) \x - y\ 



2 rT l^/'+N-l 



\Xt - Yt\ 



2T2 



ait r\Xt-Yt)\' 
\Xt-Yt\^ 



dt 



{PTf{y)Y = {mfixrw < {PTnx))m^'^^-^^r-\ 
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Since by (Al) and {A2) it is easy to estimate moments of R, the desired Harnack inequality 
follows immediately. 

In general, if a{t,x) depends on x, then the process Xf — Yt contains a non-trivial 
martingale term, which can not be dominated by and bounded drift. So, in this case, any 
additional bounded drift put in the equation for Yt is not enough to make the coupling 
successful before a fixed time T. This is the main difficulty to establish the Harnack 
inequality for diffusion semigroups with non-constant diffusion coefficient. 

In this paper, under assumptions {Al) and {A2), we are able to constructed a coupling 
with a drift which is unbounded around a fixed time T, such that the coupling is successful 
before T. In this case the corresponding exponential martingale has finite entropy such 
that the log-Harnack inequality holds; if further more (^3) holds then the exponential 
martingale is L^-integrable for some p > 1 such that the Harnack inequality with power 
holds. More precisely, we have the following result. 

Theorem 1.1. Let cr(t, x) and b(t,x) either be deterministic and independent oft, or 
satisfy (AA). 



;i) IfiAl) and {A2) hold then 



Kt\x - y\ 



Prlog/(y) < \ogPTf{x) + ^y2^l^_^.KrTy / > 1, ^, 1/ ^ T > 0. 

(2) //(Al), {A2) and{A?,) hold, then for p> {l + ^f anc? (5p,T := max{(5T, ^(v^-1)}, 
the Harnack inequality 

(TD f( \\p ^ (TD fP( w KTy/p{y/p-l)\x-y\'^ 
{PTf{y)r < {PTf^{x)) exp 



holds for all T > 0,x,y eR'^ and f e ^+{R'^). 

Theorem 11.1( 1) generalizes a recent result in |20] on the log-Harnack inequality by 
using the gradient estimate on Pf. 

Let pt{x,y) be the density of Pt w.r.t. a Radon measure fi. Then according to [261 
Proposition 2.4], the above log-Harnack inequality and Harnack inequality are equivalent 
to the following heat kernel inequalities respectively: 
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PT{y,z) 



2X^1 - e-^T 



K\x — y 



G R^T > 



and 



(1.4) 




U6MVP + 1) [( - 1) At - 6,,t] (1 - e"^-^) J ' 



x,y G M^T > 0. 



So, the following is a direct consequence of Theorem 11.11 

Corollary 1.2. Let a(t,x) and b(t,x) either be deterministic and independent oft, or 
satisfy {AA). Let Pt have a strictly positive density pt{x,y) w.r.t. a Radon measure fi. 
Then (Al) and {A2) imply (O), while {A1)-{A3) imply ( [77D . 

Next, by standard applications of the Harnack inequality with power, we have the 
following consequence of Theorem 11.11 on contractivity properties of Pt- 

Corollary 1.3. Let a{t,x) and b{t,x) be deterministic and independent oft, such that 
{A1)-{A3) hold for constant K, X and 6. Let Pt have an invariant probability measure fi. 

(1) If there exists r > K~^/X^ such that /i(e''' '^) < oo, then Pt is hypercontractive, i.e. 
||-Pt||L2(/x)~i>L4(/^) = 1 holds for some t > 0. 

(2) //|U(e^l'l^) < oo holds for allr > 0, then Pt is supercontractive, i.e. H-P* ||l2(^)^l4(^) < 
oo holds for all t > 0. 

(3) IfPte'^^'^^ is bounded for any t,r > 0, then Pt is ultracontractive, i.e. \\Pt\\L^{p)^L°°iti) < 
oo for any t > 0. 

Remark 1.1. To see that results in Corollary 11.31 are sharp, let Pt be symmetric w.r.t. 
fi. Then the hypercontractivity is equivalent to the validity of the log-Sobolev inequality 



/i(f log/2) < c/i(r(/,/)), / G cr(M'^),/i(/2) = 1 



for some constant C > 0. Moreover, if there exists a constant R > such that 
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(1.5) r(/,/)<i?2|v/|2, feC^iR'), 

we have pa{x,y) > R^^\x — y\. So, by the concentration of measure for the log-Sobolev 
inequahty, the hypercontractivity imphes /i(e^'''^) < oo for some r > 0, while the super- 
contractivity imphes /i(e^'''^) < oo for ah r > 0. Combining this with CoroUary II. 3^ we 
have the foUowing assertions under conditions {A1)-{A3) and fll.Sp : 

(i) Let K < 0. Then Pt is hypercontractive if and only if /^(e''' '^) < oo holds for some 
r > 0; 

(ii) Pt is supercontractive if and only if /i(e^l l^) < oo holds for all r > 0; 

(iii) Pt is ultracontractive if and only if P^e''' '^ is bounded for any t,r > 0. 

Therefore, conditions in Corollary 11.3( 2) and (3) are sharp for the supercontractivity and 
ultracontractivity of Pt- Moreover, as shown in that when a is constant, the sufficient 
condition /i(e''' '^) < oo for some r > /\^ is optimal for the hypercontractivity of Pt- 
So, Corollary 11.3( 1) also provides a sharp sufficient condition for the hypercontractivity 
olPf 

We will prove Theorem 11.11 and Corollary 11.31 in the next section. In Section 3 we 
extend these results to SDEs on Riemannian manifolds possibly with a convex boundary. 
Finally, combining results in Section 3 with a conformal change method introduced in [21] , 
we are able to establish Harnack inequalities in Section 4 for the Neumann semigroup on 
a class of non-convex manifolds. 



2 Proofs of Theorem 11.11 and Corollary 11.3 

Let x,y ^ M.^, T > and p > (1 + St/XtY be fixed such that x ^ y. We have 



For e e (0,2), let 



6 = V^(i-e^-(*-^)), te[o,r]. 

Then ^ is smooth and strictly positive on [0, T) such that 



(2.2) 2-KTit + i't = e, tG[0,T]. 

Consider the coupling 

dXt = a{t, Xt)dBt + h{t, Xi)dt, Xo = x, 
^^■^^ dFi = a(t, Yt)dBt + 6(t, Fi)dt + ^a(t, Fi)a(t, Xj-^(Xi - Fjdt, = y. 

Since the additional drift term ^^^a{t^y)a{t^x)~^{x — y) is locally Lipschitzian in y if 
(/14) holds, and continuous in y when a and 6 are deterministic and time independent, 
the coupling (Xf, Yt) is well defined continuous process for t < T A where C is the 
explosion time of Yt] namely, C = lining ooCn for 

Cn :=inf{t G [0,T) : \Yt\>n}, 

where we set inf = T. Let 



If C = T and 



Rs := exp 



ABt = dBt + ^a(t, XO-'(Xt - Ft)dt, t < T A C. 



is a uniformly integrable martingale for s G [0,T), then by the martingale convergence 
theorem, Rt := limtf-T-R* exists and {Rt}t&[o,T] is a martingale. In this case, by the 
Girsanov theorem {Bt\te[o,T) is a ci-dimensional Brownian motion under the probability 
Rt^. Rewrite (Q as 

, , dXf = ait, Xt)dBt + b(^, Xt)dt - ~ ^M t, Xo = a;, 

(2.4) 6 

dYt = a{t,Yt)dBt + h{t,Yt)dt, Y^ = y. 
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Since ^^^dt = oo, we will see that the additional drift — ^*^^* dt is strong enough to 
force the coupling to be successful up to time T. So, we first prove the uniform integrability 
of {-RsAc}se[o,r) w.r.t. P so that Rtac '■= ^^^s^t Rsac exists, then prove that C = T Q-a.s. 
for Q := Rt AC^ so that Q = Rt^. 
Let 

r„ = inf{tG [0,T): \Xt\ + \Yt\ >n}. 
Since Xt is non-explosive as assumed, we have r„ t C n f oo. 

Lemma 2.1. Assume (Al) and {A2). Let 9 e (0, 2),x, y G R'^ and T > be fixed. 

(1) There holds 



sup ERsAm log RsAm < 



smnn>i " " " - 2A^e(2 -e){i- e-^TT^\- 

Consequently, 

RsAC ■= hm _RsAT„A(T-i/n), s G [0,T], -Rtac •= lim-RsAC 
exist such that {RsA(}se[o,T] is a uniformly integrable martingale. 
(2) Let Q = Rtac^- Then Q(C = T) = 1 so that Q = RtF. 
Proof. (1) Let s G [0, T) be fixed. By (El), (^1) and the Ito formula, 
d\\Xt -Ytf< 2{{a{t, Xt) - ait, Yt)){Xt - Y^), dB^) + KrlX^ - Y^l'dt - ^\X, - Y^l^t 



holds for t < s A r„. Combining this with (12. 2 p we obtain 



d^ 



\Xt-Yu 



|2 



(2.5) < ^{{a{t,Xt) - cx{t,Yt)){Xt - Yt),dBt) - (2 - Kr^t + Qdt 



2 



= h{a{t,Xt) - a{t,Yt)){Xt - Yt),dBt) - ^\Xt - Yt\^dt, t<sATn. 
Multiplying by 4 and integrating from to s A r^i, we obtain 
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Jo "^t Jo '^^t uc,t '^qo 

By the Girsanov theorem, {Bt}t<T„As is the li- dimensional Brownian motion under the 
probabihty measure RsArn^- So, taking expectation Es^„ with respect to Rsatu^, we 
arrive at 

?2 dt< — — 

4t t^^o 



(2.6) E,,„ / ' d^< , se[0,T),n>l. 



By (A2) and the definitions of Rt and Bt, we have 



Jo St Jo St 

2 



<-/ i((7(t,X,)-^(X,-F,),d^,) + -^ r dt, r<sAr, 

St ^^T Jo St 

Since {i?t} is the d- dimensional Brownian motion under Rsat„^ up to s A r„, combining 
this with fl2.6D we obtain 



\x-y 



2 



ER,,At„ log RsATn = ^s,n log RsATn < ryx2nc ' S e [0, T), n > 1. 

By the martingale convergence theorem and the Fatou lemma, {-R^ac • s G [0,T]} is a 
well-defined martingale with 

Ei?.., log i?..c < 1^ = ^^To0^0^r^y s G [0, T]. 

To see that {Rsa( : s G [0,T]} is a martingale, let < s < t < T. By the dominated 
convergence theorem and the martingale property of {Rsatu ■ s G [0,T)}, we have 



E{RtAc\'^s) =E( lim i?tAr„A(r-l/n) l^s ) = hm EiRtATnA{T~l/n)\'^s) 
= lim RsATn = RsAC- 
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(2) Let (T„ = inf{t > : \Xt\ > n}. We have cr„ t P-a.s and hence, also Q-a.s. Since 
{Bt} is a Q-Brownian motion up to T A (, it follows from (12 .Sp that 



Xt AamAcS ^ \x-y\ 



6 



holds for all n > m > and t G [0, T). By letting first n f oo then m f oo, we obtain 
Q(C < t) = for all t e [0,T). This is equivalent to Q(C = T) = 1 according to the 
definition of C. □ 

Lemma [2m ensures that under Q := Rtai^^, {Bt}te[o,T] is a Brownian motion. Then by 
(12. 4p . the coupling {Xt, Yt) is well-constructed under Q for t G [0, T]. Since ^^^dt = oo, 
we shall see that the coupling is successful up to time T, so that Xt = Yr holds Q-a.s. 
(see the proof of Theorem 11.11 below). This will provide the desired Harnack inequality 
for Pf as explained in Section 1 as soon as Rta<; has finite p/{p — l)-moment. The next 
lemma provides an explicit upper bound on moments of Rta(- 

Lemma 2.2. Assume {A1)-{A3). Let Rt and he fixed for 9 = Oj-- We have 



(2.7) sup E<^ i?s/\<^ exp 
se[o,T] 

Consequently, 



6'^ 

85^ Jo 



< 



exp 



9tKt\x — 



L452(2-0r)(l-e-^^^) 



(2.8) 

holds for 



sup ^R]Xr < exp 



se[o,T] 



8(52 (2 _ ^ eT\T){l - e-^TT^ 



46^ + AOtXt^t ' 



Proof. Let 6 = 6t- By (12. 5p . for any r > we have 



^s,n exp 

< exp 

< exp 





r\x — y\ 



dt 



-E^^^exp 
tKtIx — 1/|2 



2r 



1 



Ler(2-er)(l-e-^^^ 



Es,n exp 



((a(t,X0-a(t,r,))(X,-F,),d5,) 

1/2 



fl2 



dt 
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where the last step is due to {A3) and the fact that 



for a continuous exponential integrable martingale Mf. Taking r = 6^/(86^), we arrive at 



E^.n exp 



9^ r'^'-ix, 



8^T Jo 



-dt 



< 



6tKt\x — y\ 



L452(2-^r)(l-e-^^'r)J 



n > 1. 



This implies (12 .yp by letting n — )■ oo. 

Next, by (^42) and the definition of we have 



(2.9) 



Es.n exp 



Jo Si 

\ait,Xty'iXt-Yt)\' 



+ 



2 Jo 



dt 



Noting that for any exponential integrable martingale Mf w.r.t. Rsf^rJ^i one has 



E,,„exp[rTMt + rT(M)t/2] = E,,„exp[rTM4 - rlq{M)t/2 + rT(grT + l){M)t/2] 
< (E,,„exp[rrgMt -4g2(M)i/2])'/''(E,,„exp 



rrg(rrg + 1) ^^^^^ J \ (9-1)/^ 



exp 



TTqirxq + 1) 



2(g-l) 
it follows from ([23]) that 



g > 1, 



(2.10) 
Take 



^RlX:: < E.„exp 



2(g-l)A 



jAt„ 



\Xt-Y,\^ 



dt 



(2.11) 



g = 1 + J 1 + 1 



which minimizes q{qrT + l)/(g — 1) such that 
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(2.12) 



qrriqrT + 1 + v^ry + 1 . ^- — . 

ox2^^ IN = / + 1 + VrT(rT + 1)) 



2 Aji 



+ 
86^' 



Combining f l2.10p with f l2.7p and f l2.12p . and noting that due to fl2.1ip and the definition 
of Tt 



we obtain 



1 + r, 



-1 

T 



26r + OtX 



25x + 29j'Xt 



Ei?^+;j < exp 



-852(2 - ^r)(5r + ^tAt)(1 - e-^^^) 
According to the Fatou lemma, the proof is then completed by letting n ^ oo. 



□ 



Proof of Theorem \l.l[ Since (^43) also holds for Sp^T in place of St, it suffices to prove the 
desired Harnack inequality for 6t in place of 6p^T- 



(1) By Lemma [2?H {-RsAc}se[o,T] is an uniformly integrable martingale and {Bt}t<T is 
a (i-dimensional Brownian motion under the probability Q. Thus, Yt can be solved up to 
time T. Let 

r = inf{tG [0,T] : Xt = Yt} 

and set inf = oo by convention. We claim that r < T and thus, Xt = Yt, Q-a.s. 
Indeed, if for some w G f2 such that t{uj) > T, by the continuity of the processes we have 



inf \Xt-Yt\\u)>0. 

ie[o,T] 



So, 



T IV V12 

dt = OO 



\Xt~YA 
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holds on the set {r > T}. But according to Lemma [2.21 we have 



dt < oo, 



'0 

we conclude that Q(r > T) = 0. Therefore, Xt = Yx Q-a.s. 

Now, combining Lemma 12.11 with Xj- = It and using the Young inequality, for / > 1 
we have 



Prlogfiy) = EQ[log/(FT)] = E[i?TAC log /(X^)] 

< Ei?TAclog^TAC + logRfiXr) < log Prfix) + 



Kt\x - y\ 



2Xle{2-e){l-e~^TTy 

This completes the proof of (1) by taking 6 = 1. 

(2) Let 9 = 9t- Since Xt = and {i?(}(g[o,T] is the ci- dimensional Brownian motion 
under O, we have 



(2.13) {PrfiyW = mfiYrW = (E[i?rAc/(^r)])^ < {PTFixmERf^-'^^^-' 



^TAC 



Due to (12.11) we see that 



p — 1 ' A5t{5t + OtXt) 
So, it follows from Lemma 12.21 and (12. ip that 



(EH: 



p/(p-i)np-i 

TAC 



(Ei?^+"^^)f-^ < exp 



exp 



(p - 1)^t/'S^t(25t + OtXt) \x - 
852(2 - eT){6T + ^tAt)(1 - e-^TT^ 
KT^/piVP-l)\x-y\^ 



AST[{Vp-l)^T-ST]{l-e-KTT) 
Then the proof is finished by combining this with ( 12.13^ . 



□ 



Proof of CorollarylTE Let / G ^^{R'^) be such that //(/p) < 1. Let p > (1 + 5/\f. By 
Theorem 11.1( 2). we have 
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where 5p = max{5, f !)}■ Integrating w.r.t. /i(dx) and noting that /i is Pj-invariant, 

we obtain 



(2.14) 



iPtf{y)r / exp 



K^i^-l)\x-y\ 



fi{dx) < 1. 



45[(^- 1)A-(5](1 -e-^*)- 

Taking f = n A {pt{y, ■)Y^^ and letting n f oo, we prove the first assertion. 

Next, let 5(0,1) = {x G M'' : \x\ < 1}. Since fi is an invariant measure, it has a 
strictly positive density w.r.t. the Lebesgue measure so that fi{B{0, 1)) > (cf. [6]). Let 
p > (1 + 26/ \y. We have Sp = - l)A/2 and thus. 



46p[{^-l)\-6p 
Combining this with f l2.14p and noting that 



A2(^-l) 



K^{^-l){l + \y\y 



> /i(E(0, l))exp 



4(5[(v^- 1)A-5](1 -e-^*)J 



we obtain 



(2.15) 



iPtfiy)r<C,exp 



K^{l + \y\y 



t>0,ye 



for some constant Ci > and all / G ^^{W'-) with n{fP) < 1. Since 



lim lim 



t^Z A2(^ - 1) (1 - e--^*) A2 ' 
for any r > K'^/X'^ there exist p > {1 + 26t/\Y, > 1 and ti > such that 



{PtJ{y))^''<C2e^^'^\ yeR'Je 



r\y\^ 
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holds for some constant C2 > 0. Thus, /i(e''' '^) < 00 imphes that \\Pti\\LP{fi)^LPi^{/i) < 00. 
Since ||Ps||l'j{^j) = 1 holds for any q G [l,oo], by the interpolation theorem and the 
semigroup property one may find ^2 > ti such that 

(2-16) ||-Pt2l|L2(M)^L4(^) < 00. 

Moreover, by [131 Theorem 3.6(ii)], there exist some constants 77, C3 > such that 

\\Pt- f^\\LH^.)<C3e-'^\ t>0. 

Combining this with fl2.16p we conclude that H-P* ||l2(^)->l4(/^) < 1 holds for sufficiently 
large t > 0, i.e. (2) holds. 

Finally, (3) and (4) follow immediately from f l2.15p and the interpolation theorem. □ 

3 Extension to manifolds with convex boundary 

Let M be a rf- dimensional complete, connected Riemannian manifold, possibly with a 
convex boundary dM. Let N be the inward unit normal vector filed of dM when dM ^ 0. 
Let Pt be the (Neumann) semigroup generated by 

L ■.= ij\A + Z) 

on M, where ip G C^(M) and Z is a vector field on M. Assume that ip is bounded 
and 

(3.1) Ric-VZ >-Ko 

holds for some constant Kq > 0. Then the (refiecting) diffusion process generated by L is 
non-explosive. 

To formulate Pt as the semigroup associated to a SDE like (11. ip . we set 



(3.2) 
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Let d/ denote the Ito differential on M. In local coordinates the Ito differential for a 
continuous semi- martingale Xt on M is given by (see [10] or [1]) 

1 

(d,x,)' = dx^ + rf,(x,)d(x\x^),, i<k<d. 

Then is the semigroup for the solution to the SDE 



(3.3) diX, = a(Xi)$idSi + h{X^)dt + iV(Xi)d/t, 

where Bt is the (i-dimensional Brownian motion on a complete filtered probability space 
(f2, {t^t}t>o,P), $t is the horizontal lift of Xt onto the frame bundle 0(M), and /t is the 
local time of Xt on DM. When DM = 0, we simply set k = 0. 

To derive the Harnack inequality as in Section 2, we assume that 



(3.4) A := inf a > 0, 6 := sup a — inf a < oo. 

Now, let x,y & M and T > be fixed. Let p be the Riemannian distance on M, i.e. 
p{x, y) is the length of the minimal geodesic on M linking x and y, which exits if dM is 
either convex or empty. 

Let Xt solve (13. 3p with Xq = x. Next any strictly positive function ^ G C([0,T)), let 
Yt solve 



djYt = a{Yt)Px„Y.^tdBt + b{Xt)dt - ^i^iM^i^Vp(X„ ■){Yt)dt + N{Yt)dk 

for Yq = y, where k is the local time of Yt on dM, and PxtXt '■ '^XtM Ty^M is the 
parallel displacement along the minimal geodesic from Xt to Yt, which exists since dM is 
convex or empty. As explained in [H Section 3] , we may and do assume that the cut-locus 
of M is empty such that the parallel displacement is smooth. Let 

dBt = dBt + Vp(-, Yt){Xt)dt, t < T. 
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By the Girsanov theorem, for any s G (0, T) the process {Bt}t£[o,s] is the d-dimensional 
Brownian motion under the weighted probabihty measure -R^P, where 



I III /■X I 111 11(11/1111^^ 

(3.5) Rs := exp 

Thus, by (13. 2 p we have 



djXt = V2i,{Xt)<^tdBt + {ij'Z)iXt)dt - ^t^l^Vp(-, Yt)iXt)dt + N{Xt)dk, 

lit 

djYt = V2^{Yt)^tdBt + {iP^Z){Yt)dt + N{Yt)dk. 
Let C, G C^([0,T)) be strictly positive and take 

Repeating the proof of (4.10) in [27] we obtain 



dp{X„Y,) < (a(X,)-a(F0)(Vp(-,Ft)(Xi),$id5,) + Kip(X„Fi)dt-^^^^dt, t<T, 
where 



^i = ^o||V^|lL + 2||^||oo||V^|| 

This imphes that 



^p{X^ < |-p(X„r,)(a(X,)-a(F,))(Vp(-,r*)(X,),$td5,)-^^%^(2-J^6+eOdt 
holds for t < T and 



K:=2Ki + ||Va||L 

3.6 

= 2Ko\ml + 4||Z|U||VV'||oo||^||oo + 2||VVir 



ooWY \\oo- 



oo" 



In particular, letting 

18 



6 = ^(i-e^(*-^)), te[0,T],9e{0,2), 



we have 



2-K^, + C = e. 

Therefore, the following result follows immediately by repeating calculations in Section 2. 
Theorem 3.1. Assume that dM is either empty or convex. Let (\4-i^ and Z, be bounded 



such that 

:= 2i^o||^|lL + 4||^||oo|| V^IUII^Iloo + 2|| VV^IIL < oo. 

Then all assertions in Theorem and Corollaries hold for Pt the {Neumann) 

semigroup generated by L = tp'^{A + Z) on M , and for constant functions K,6 := sup ip — 
iniip and X := inf 

4 Neumann semigroup on non-convex manifolds 



Following the line of [21], we are able to make the boundary from non-convex to convex 
by using a conformal change of metric. This will enable us to extend our results to the 
Neumann semigroup on a class of non-convex manifolds. 

Let dM 7^ with the inward normal unit vector field. Then the second fundamental 
form of dM is a two-tensor on the tangent space of dM defined by 

I(X, Y) := -{VxN, Y), X,Y e TOM. 
Assume that there exists k > and Kq G M such that 

(4.1) mc-VZ>-Ko, I>-K 

holds for M and a vector field Z. We shall consider the Harnack inequality for the 
Neumann semigroup Pt generated by 
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L = A + Z. 

To make the boundary convex, let / G C^(M) such that / > 1 and A^log/|aAf > '^• 
By Lemma 2.1], dM is convex under the metric 

(■,-)' = r^(-,-)- 

Let A' and V be the Laplacian and gradient induced by the new metric. We have (see 
(2.2) in PI]) 

L = f-\A' + Z'), Z' = fZ + ^Vf. 

Let Ric' be the Ricci curvature induced by the metric (■, ■)'. We have (see the proof of 
[271 Theorem 5.1]) 

Ric'- V'Z' > -Kf {■,■)' 

for 

(4.2) Kf = sup {Kf - dAf + id- 3)|V/p + 3|Z|/|V/|}. 

Applying Theorem 13. II to the convex manifold (M, (-, ■)'), ip = and 

K = 2Ki\\f-'\\^ + 4||z'ir^ii v'ri'ooiirloo + 2|v'ri'oo' 

< 2Kl + 4\\fZ + {d- 2)V|UII V/IU + 2|| V/IIL, 

where || ■ ||' is the norm induced by (■, ■)' and we have used that / > 1, we obtain the 
following result. 

Theorem 4.1. Let ( |^. ip hold for some k > and Kq G M, and let Pt he the Neumann 
semigroup generated by L = A + Z on M. Then for any f G C^{M) such that inf f = 1, 
N log f\dM > and K < oo, where K is fixed by l \4-S^ and all assertions in Theorem 

\1.1\ and Corollaries \1.2\ and \1.3\ hold for constant functions K, 6 := sup — inf f~^, and 
A :=inf/-i. 
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Remark 4.1 A simple choice of / in Theorem 14.11 is / = o p^, where pQ is the 
Riemannian distance to the boundary which is smooth on {pQ < r^} for some r^^ > 
provided the injectivity radius of the boundary is positive, and / G C^([0, oo)) is such 
that /(O) = l,/'(0) = K and /(r) = /(r^-) for r > r^. In general, / is taken according 
to tt and bounds of the second fundamental form and sectional curvatures, see e.g. 
pl| [27] for details. With specific choices of /, Theorem 14.11 provides explicit Harnack 
type inequalities, heat kernels estimates and criteria on contractivity properties for the 
Neumann semigroup on manifolds with non-convex boundary. 

Acknowledgement. The author would like to thank the referees and Dr. Wei Liu for 
useful comments and corrections. 
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